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ABSTRACT 



We show that the radio, X-ray and 7-ray spectrum of the supernova remnant 
RX J1713. 7-3946 can be accounted for with the simplest emission model, where 
all of these emissions are attributed to a population of relativistic electrons inter- 
acting with the cosmic microwave background radiation, IR interst ellar photons, 
and a background magnetic field. With a spectral inversion method (IJohns fc Lin 
19921 ). the parent electron distribution and its uncertainties are derived from the 
observed photon spectrum. These results are independent on the model of parti- 
cle acceleration and strongly support the leptonic scenario for the TeV emission. 



Subject headings: Acceleration of particles — ISM: supernova remnants 
ods: miscellaneous — Radiation mechanisms: non-thermal 



Meth- 



1. INTRODUCTION 

Relativistic cosmic rays (CRs) with energies lower than the location of the "knee" of 
the cosmic ray spectrum (~ 10 15 eV) are com monly believ ed to be accelerated at the shock of 
supernova remnants (SNRs) in our Galaxy (jHillas I l2005h . Evidence of particle acceleration 
by shocks of SNRs first comes from radio observations of the remnants, where the radio 
emission is produced by relativistic electrons via the synchrotron process. The discovery of 
synchro tron X-ray emission f rom SNR 1006 reveals acceleration of TeV electrons by these 
shocks ( iKoyama et al. Ill995l ). However, direct evidence of proton and ion acceleration up 
to the spectral "knee" remains elusive. 
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rays from SNRs can be produced via a hadronic process, where neutral pions produced by 
inelastic collisions of energetic hadrons decay into 7-rays . Energetic leptons can produce 
7-rays through inverse Compton (IC) scattering of low frequency background photons. Both 
energetic leptons and hadro ns can also produce 7-rays through the bremss trahlune process 



(IBlumenthal fc Gouldlll970h . which is negligible for SNR RX J1713. 7-3946 (lAharonian et al. 
20061 ). Multi- wavelength observations of photon emission from this SNR in combination 



with theoretical considerations and/or detailed nu merical modelings have been used to areue 
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explorations, the dominant TeV emission mechanism remains a matter of debate. Future 
neutrino experiments are expected to measur e contributions to the observed TeV emission 



via the hadronic process (jVissani et al. 1 1201 if ). 



SNR RX J1713. 7-3946 has been extensively studied from radio, IR (Acero et al. 2009), 
X-ray (Koyama et al. 1997; Uchiyama et al. 2003; Cassam-Chena et al. 2004), GeV 7- 
ray (Abdo et al. 2011), to TeV 7-ray band (Muraishi et al. 2000; Enomoto et al. 2002; 
Aharonian et al. 2006). The TeV 7-ray observations reveal a shell- like structure matching 
closely to the radio and non-thermal X-ray shells (Aharonian et al. 2006). Most recently, 
observations made with the Fermi space telescope show that the spectrum of this source in 
the GeV band is very hard with a power-law photon index of T = 1.5 ±0.1, which is difficult 
to accommodate in a hadronic scen ario, but agrees well with the IC origin of 7-rays in the 
leptonic scenario (lAbdo et al.ll201ll ). On the other hand, b y considering the po tentially high 
inhomogeneity of the shocked interstellar medium (ISM), llnoue et al. I (120111 ) argue that a 
hadronic model may still explain the Fermi observation. 

However, even if future observations reveal 7-rays via the hadronic process from isolated 
SNRs, this does n ot imply that they dominate the Galactic CR flux observed near the Earth. 
Ave et al. I (120091 ) recently showed that the de-propagated Galactic CR source spectrum is 
much softer than that given by most diffusive shock acceleration models. Furthermore, it 
has already been noticed that most of the SNRs are overlapping in superbubbles, where 
the partic le acceleration mechanism may be quite different from that operatine in isolate d 
remnants JHiedon et alill998l : IParizot et al. Ibood ; butt H2OO9I : iFerrand & Marcowith Iboioh . 



Given the good spatial correlation between images made at different energy bands (Acero 
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et al. 2009), certain distributions of energetic particles or mechanisms of particle accelera- 
tion are usually introduced to fit the spatially integrated broadband spectrum in all these 
previous studies. Different scenarios of the background radiation field and ISM are also con- 
sidered. Uncertainties in the particle acceleration process at the SNR shocks, especially in the 
spatial diffusion coefficient of energetic particles and in jection mechanis ms of supra-thermal 
particles, have made the relevant studies inconclusive ( )Fan et al. II2010I ). 

In this Letter, we show that the current observations, though very extensive, may not 
justify these sophisticated modelings. Assuming that the TeV 7-rays originate from the 
IC by relativistic electrons of the cosmic microwave background radiation (CMB) and an 
IR background with a black-body spectrum suitable for the Galactic environment of the 
remnant, we show that one can derive the parent electron distribution and its uncertainties 
directly from 7-ray observations with a well-established spectral inversion method (§ 
It is also shown that by adjusting the magnetic field and the extrapolation of the parent 
electron distribution toward low- and high-energies, the radio to X-ray spectrum can also 
be reproduced. Therefore the simplest emission model, where the broadband emission of 
SNR RX J1713. 7-3946 is attributed to a population of relativistic electrons interacting with 
a background magnetic field, the CMB and an IR photon background, can fully account for 
the radiation spectrum and gives direct constraints on the particle acceleration processes 
(§ E]). Although these results do not preclude some level of proton 7-rays being masked by 
the more intense electron 7-rays , improved spectral measurements are needed to go beyond 
such a simple model. Conclusions and implications of our study are drawn in § HJ 



2. SPECTRAL INVERSION METHOD 

The primary goal of the study of astrophysical sources is to use characteristics of the 
observed emissions to probe the underlying physical processes. In general, the related ra- 
diation mechanisms are the first and most important processes to be studied. When the 
quality of the observational data is not sufficiently good, some forward modeling approach is 
usually taken to explain the relevant observations. The consequent results are usually model 
dependent since different assumptions can be introduced in the proposed models. 

When the relevant radiative mechanisms are well established and the observed data 
have a high quality, a spectral inversion method may be used to derive the distribution of 
the parent particle distribution from the observed radiation spectrum. The parent particle 
distribution is more directly connected to the underlying physical processes than the observed 
radiation spectru m. This is a classi cal inversion problem and has been studied extensively 



in solar physics ( IJohns fc Linlll992l ). In the hadronic scenario for the TeV emission from 
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shell-type SNRs, IVissani I ( 120061 ) shows that the observed TeV emission spectrum can be 
used to derive the neutrino spectrum from these sources. However, it is well known that the 
leptonic scenario for the TeV emissi on from SNRs has fewer model parameter s and faces less 



challe nges than the hadronic one (ILiu et al. I l2008t iKatz fc Waxman 1 120081 ; Ellison et al. 



2010|). In light of recent simplifications of the treatment of the IC process (IPetruk 



2009 



the spectral inversion method can be readily applied to the leptonic scenario for the TeV 
emission from a few shell type SNRs. In the following, we will only consider the CMB and 
an IR background as seed photons for the IC process. 



2.1. IC Radiation of Mono-energetic Electrons in a Black-body Photon Field 



Petruk I (120091 ) has explored the IC radiation of mono-energetic electrons in black-body 
photon fields. For the sake of completeness, we summarize the key results relevant to our 
study in the following. The spectral distribution of the volume emissivity of an isotropically 
distributed electron population due to the IC process is given by (Blumenthal & Gould 1970) 

P(k) = ck f dyN{y) / den ph {e)a lc {k,e;y) (1) 



where e and k are the photon energies before and after the interaction, respectively, c, 7, 
N(y), n p h(e), and a\c are the speed of light, the Lorentz factor, spectral distribution of 
electrons, the energy distribution of the background photons, and the angle-integrated IC 
cross-section, respectively. Equation ^ can be rewritten as 

P(k) = J dyN(y)p(y,k), (2) 

where ^(7, k) represents the spectral distribution of IC radiation of mono-energetic electrons 
with a Lorenz factor 7. 

In some astrophysical cases, n p h(e) can be represented by the isotropic blackbody radia- 
tion n p h(e) = 87re 2 //i 3 c 3 [exp(e/e c ) — 1] where e c = k B T, and k B , h, and T are the Boltzmann 
constant, the Planck constant, and the temperature of the radiation field, respectively. Then 
p(y, k) in Equation (J2]) can be approximated as 

P(7,t)^fc») (3) 
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where r] c = e c k / (m c c 2 ) 2 and 770 = k 2 / '[A^m e c 2 ''( , ym e c 2 — k)}, m e is the electron rest mass 
(Petruk 2009). With the approximation above, the double integral in Equation ([1]) for 
calculating the IC emissivity is simplified into a single integral as given by Equation (J2J). 
This simplification facilitates further steps in our spectral inversion method significantly 
When there are multi-blackbody components, one can obtain an I for each temperature and 
add them together with the appropriate weight of their photon energy density to get the 
overall p. 



To derive the electron distribution ^(7) from the observed photon spectrum, which 
is proportional to P(k), one needs to use energy bins of the observed photon spectrum to 
discretize Equation (Johns & Lin 1992): 



where k n corresponds to the highest energy of the observed photons, 7, is the minimum 
Lorentz factor of electrons that can contribute to the IC 7-rays with an energy 
which is what we will obtain directly from the inversion method, is the electron density in 
the Lorentz factor interval (jj, jj+x) averaged over p(j, ki). 

So, for the n energy bins of the observed photons, P(ki) forms annxl matrix, and the 
de-convolved electron distribution, Nfy), forms another nxl matrix. Then Equation (jSJ) 
can be expressed in a matrix form: 



2.2. 



Matrix Formulation of the IC Emission 



p ( k i) = E ( / # (7)P(7, ki)(h) = E / P(7, fci)d7) (5) 





(6) 



where the IC power matrix is given by: 




(7) 



Then the electron distribution can be derived by inverting this matrix equation 



n 
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^^In practice, at a given Lorentz factor 7, p(fc) may peak at several photon energies if there are multi- 
blackbody components. We adjust 7 so that the highest photon energy, where p(k) peaks, is equal to fej. 
Then 7^ is chosen as a factor of 1.3 lower than this Lorentz factor. 
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where the matrix is the inversion of the matrix (3. 
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for j = i , 
for j < i . 
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This kind of inversion method was first introduced by iJohns fc Lid (119921 ) to derive the 
parent electron distribution from the optically thin spectrum of bremsstrahlung photons. 
The IC of blackbody spectra has similar characteristics as the bremsstrahlung process after 
adopting the approximation discussed in § 12.11 One therefore can readily adopt the formulas 
in this method. Note that is gi ven explicitly here, and there should be a minus sign in 
the original Equation (14) given by I Johns &: Lin 



Since there are uncertainties in the observed 7-ray spectrum, the propagation of errors 
from the photon spectrum to the derived electron spectrum should be considered carefully. 
Fortunately, the matrix Equation (jSJ) offers a very straightforward way of calculating the 
errors in the electron spectrum. The densities of electrons simply linearly depend on the 
photon fluxes, the errors of the electron densities are then given by: 



1/2 



(10) 



3=1 



Thus, we can use this matrix formulation to obtain the electron distribution and errors from 
the observed photon fluxes and errors. 



RESULTS 



We now apply our inversion met hod to the GeV and TeV observations of S NR RX J1713- 
3946 by the Fermi-LAT and HESS (lAbdo et al.ll201ll ; lAharonian et al.l 120061 ). respectively. 
Besides the CMB, we as sume IR photons wit h a temperature of 30 K and an energy density 



of 1 eV cm -3 following iPorter et al. I ( 120061 ). Figure [T] shows the 7-ray spectrum of this 
SNR along with the resulting electron spectrum after applying the inversion method. Due 
to large errors and fluctuations of the observed 7-ray fluxes, negative electron densities are 
inferred for some energy intervals, which is clearly not physical but mathematically expected. 
To address this issue, a running smooth of the observed fluxes is done before applying the 
inversion method and the corresponding data are indicated by the squares in the Figure. 
Even with such a smoothing, the error of the derived electron density is still large in some 
energy bins. Future observations with improved 7-ray flux density data will give a better 
electron spectral measurement. 
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To verify the validity of this method, one may calculate the IC and synchrotron spectra of 
the obtained electron distribution and compare them with observations. To do this properly, 
one first needs to interpolate the obtained electron spectrum since the energy bins are wide 
and the radiation spectra are sensitive to spectral details. We use a 2nd-order polynomial 
to do the interpolation. To calculate the synchrotron spectrum, one also needs to specify 
the magnetic field and do the extrapolation of the electron distribution toward low and high 
energies. We assume power-law distributions for these extrapolations with the power-law 
indices as free parameters. By adjusting the magnetic field and the two power-law indices, 
one can fit the observed radio to TeV spectrum, as shown in Figure [2j For the best fit, 
the value of the magnetic field B = 15 /j,G and the values of the low- and high-energy end 
electron spectral indices are 2.1 and 10, respectively. This smooth electron distribution is 
shown as the red solid line in the right panel of Figure [TJ 



4. DISCUSSIONS AND CONCLUSIONS 

With a spectral inversion method, we show full consistency of the simplest leptonic 
model for the radio to TeV spectrum of SNR RX J1713. 7-3946. The model does not prescribe 
the particle acceleration process and assumes only the CMB and an IR background as the 
seed photons for the IC scatter. The current observations, though extensive, therefore may 
not justify very sophisticated modeling of this source, and there is no prominent evidence 
of emission from energetic hadrons. The errors of the derived electron distribution are large 
and the resulting electron distribution can be readily fitted with a simple analytical function: 
7 _2 exp [—(7/I.3 x 10 7 ) ' 6 ] (blue solid line in the right panel of Fig. [1]), which also fits the 
overall radiation spectrum (blue lines in Fig. [5]). Although the difference between this 
analytical function and the distribution obtained from the inter- and extrapolation of the 
derived electron densities appears to be prominent, it is not statistically significant due these 
large errors. 

These results have not taken into account the effect of optical background photons. 
Inclusion of this component, though more involving, is straightforward. Actually, even the 
adoption of blackbody spectra is not necessary. With the blackbody background photon 
spectrum, one may obtain some analytical expressions to better appreciate these results. For 
the spectral inversion method to be applicable, one just requires that, for a given Lorentz 
factor of electrons, the IC spectrum cutoff sharply at high energies so that the IC cross section 
may be represented by a triangle matrix. To estimate the effect of the optical background 
photons, we derive the electron distribution assuming the CMB alone. With a lower overall 
photon energy density, the inferred magnetic field is only 8.5/iG, which is expected. The 
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analytical approximation of the electron distribution is very similar except that the cutoff 
Lorentz factor changes from 1.3 x 10 7 to 1.6 x 10 7 , which is caused by the higher energy 
of IR photons than that of the CMB. The spectral indexes of the electron distribution 
extrapolated to low and high energies are 2.2 and 10, respectively, which are essentially the 
same. Therefore, with an extra optical photon background, we expect that the magnetic 
field increases even more and the cutoff energy of the electron distribution decreases slightly. 



Fan et al. 
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may not be justified for the following reasons: 1) there are systematic errors in the observed 
X-ray and 7-ray fluxes; 2) the Suzaku X -ray spectrum is ob tained from a part of the remnant 



and re-scaled to give the overall flux es (Tanaka et al. | 2008l ). which may not be well justified 



due to the asymmetry of the source (lAcero et al.ll2009h : 3) the complex source structure also 



betrays the simple one zone emission model (jAcero et al. 2009 ). and both the magnetic field 



and the electron distribution may vary significantly across the remnant 0; 4) the background 
radiation field is not well known, which constitutes another systematic uncertainty; 5) the 
particle acceleration proc ess can be intrinsically complex without simple distributions of 
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the acceleration particles (IPetrosian fc Liu Il2004t iMewaldt et al. I l200ll : IChang et al. 1 12008 



Adriani et al. 



201ll ). Nevertheless, the method 



can be used to other sources where the emission is dominated by optically thin IC process. 
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Fig. 1. — Left: the observed (circle) and smoothed (square) GeV and TeV 7-ray fluxes for 
RX J1713. 7-3946. (Abdo et al. 2011, Aharonian et al. 2006). Right: the derived electron 
distribution with error bars. The red solid line represents the inter- and extrapolated electron 
distribution, which is used to calculate the synchrotron and IC radiation spectra in Figure 
HJ The blue solid line shows an analytical function, which can also reproduce the observed 
radiation spectrum. 
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Fig. 2. — Comparison of the observed radio (Acero et al. 2009), X-ray (Tanaka et al. 2008) 
and 7-ray fluxes with the synchrotron (solid) and IC (dashed) spectra of the derived electron 
distributions using our inversion method. The blue lines are for the analytical distribution, 
whose parameters are described in § [3j The red lines are for the inter- and extrapolated 
electron distribution, where the dotted and dot-dashed lines are for the IC of IR and CMB 
photons, respectively. 



